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Abstract. We first demonstrate that globular cluster Gl could have been formed by tidal interaction between M31 and a 
' nucleated dwarf galaxy (dE,N). Our fully self-consistent numerical simulations show that during tidal interaction between M31 

I and Gl's progenitor dE,N with Mb ~ -15 mag and its nucleus mass of ~ 10' Mq, the dark matter and the outer stellar envelope 

, of the dE,N are nearly completely stripped whereas the nucleus can survive the tidal stripping because of its initially compact 

nature. The naked nucleus (i.e., Gl) has orbital properties similar to those of its progenitor dE,N. The stripped stars form a 
, metal-poor ([Fe/H] 1) stellar halo around M31 and its structure and kinematics depend strongly on the initial orbit of Gl's 

progenitor dE,N. We suggest that the observed large projected distance of Gl from M31 (~ 40 kpc) can give some strong 

constraints on the central density of the dark matter halo of dE,N. We discuss these results in the context of substructures of 

M3rs stellar halo recently revealed by Ferguson et al. (2002). 

Key words, galaxies: halos — galaxies: individual (M31) 
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1. Introduction 

A growing number of photometric and spectroscopic observa- 
tions have suggested that Gl (=Mayall II), which is one of the 

I brightest globular clusters belonging to M31, has very unique 

' physical properties as a globular cluster (e.g., Meylan et al. 
2001). These include a possible intrinsic metallicity dispersion 

■ among its stellar population (Meylan et al. 2001), the large cen- 

[ tral velocity dispersion of ~ 25 km s ' (e.g., Djorgovski et al. 
1997), the very flattened shape with mean ellipticity of 0.2 and 
significantly high central surface brightness (Rich et al. 1996; 
Meylan et al. 2001), and the possible existence of a black hole 
with a mass of 20000 Mq (Gebhardt et al. 2002). In spite of 
its extraordinary nature, Gl is observed to be on the My-cro re- 
lation (where My and ctq are total magnitude in V-band and 
central velocity dispersion) defined not by elliptical and dwarf 
galaxies but by globular clusters, which implies that Gl looks 
like a genuine globular cluster (Meylan et al. 2001). 

One of the possible scenarios of Gl formation is that Gl 
is the surviving nucleus of an ancient nucleated dwarf galaxy 
with its outer stellar envelope almost entirely stripped by M3 1 's 
strong tidal field (Meylan et al. 1997, 2000, 2001). Such a sce- 
nario has already been suggested by Zinnecker et al. (1988) 
and Freeman (1993), and the viability of the scenario has been 
extensively discussed by many authors in terms of oj Cen for- 
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mation (e.g., Hilker & Richtler 2000; Dinescu 2002; Gnedin 
et al. 2002; Zhao 2002; Bekki & Freeman 2003; Mizutani et 
al. 2003). However, no theoretical attempts have been made so 
far to investigate (1) whether M3rs tidal filed is strong enough 
to transform a dE,N into Gl and (2) what observable evidence 
of the past destruction of Gl's progenitor dE,N we can find in 
the M31 halo regions. The above point (2) is very important, 
because Ibata et al. (2001) and Ferguson et al. (2002) have re- 
cently discovered M3 1 's stellar halo substructures, which could 
have been formed from tidal destruction of M3rs satellite 
dwarfs. 

In this paper, by using numerical simulations, we first 
demonstrate that Gl can be formed from a dE,N during tidal in- 
teraction between the dE,N and M3 1 . Our fully self-consistent 
numerical simulations demonstrate that the stellar envelope of 
dE,N with Mb ~ -15 can be nearly completely stripped by 
the strong tidal field of M3 1 whereas the central nucleus can 
remain intact owing to its compactness. We suggest that this 
naked nucleus orbiting M31 is a giant globular cluster (i.e., 
Gl). The morphological transformation from dE,Ns into very 
compact stellar systems was originally investigated by Bekki 
et al. (2001, 2003) for ultra-compact dwarfs labeled as "UCD" 
(Drinkwater et al. 2003) and called "galaxy threshing" (Bekki 
et al. 2001), though they suggested that galaxy threshing is also 
important for the formation of giant globular clusters such as a> 
Cen and Gl. The present study confirms this early suggestion 
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and discusses the physical relationship between UCDs, oj Cen, Table 1. Model parameters 
andGl. 



2. Model 

Since our numerical methods and techniques for modeling col- 
lisionless, self-gravitating systems of dE,Ns already have been 
described in detail by Bekki et al. (2003), we give only a brief 
review here. We consider that a dE,N with a mass and size sim- 
ilar to that observed for the dE,N types orbiting M3 1 . The dE,N 
is modeled as a fully self-gravitating system and is assumed to 
consist of a dark matter halo, a stellar component and a nu- 
cleus: We first investigate this "three-component" dE,N model. 
For convenience, the stellar component (i.e., the main baryonic 
component) is referred to as either the "envelope" or the "stel- 
lar envelope" so that we can distinguish this component from 
the stellar nucleus. The density profile of the dark matter halo 
with the total mass of M^m in the dE,N is represented by that 
proposed by Salucci & Burkert (2000): 



Pdm(r) = 



Pdm,Q 



(r + fldm)('-^ + adm)^ ' 



(1) 



where pdm,o and Odm are the central dark matter density and the 
core (scale) radius, respectively. This model profile is consis- 
tent with observations and different from the predictions of the 
standard CDM model (Navarro, Frenk, & White 1996, here- 
after NEW). The dark matter core parameters, pdm,o, Odm, and 
Mo (where Mq is the total dark matter mass within a^) have a 
clear observed correlation, Mq = 4.3 x 10'^(^)^^^Mo (Burkert 
1995). 

The mass (luminosity) and the scale length of the stellar en- 
velope of the dE,N is modeled according to the observed scal- 
ing relation of Ferguson & Binggeli (1994): 



logadw[pc] = -0.02Mb -I- 2.6 



(2) 



for faint dwarfs (Mb > -16), where Odw and Mb are the scale 
length of the exponential profile and the absolute B-band mag- 
nitude, respectively. The projected density of the envelope with 
Mb and the total mass of Mdw (and Mdw/^B = 2) is represented 
by an exponential profile with a scale length Udw The projected 
density profile of the nucleus with mass M„ is represented by a 
King model (King 1964) with a core radius of and a central 
concentration parameter c of 1 .0. 

The nuclei typically contribute about a few percent of the 
total light of dwarfs (Binggeli & Cameron 1991; Freeman 
1993) and the present-day mass of Gl is estimated to be 1.5 x 
10^ Mq for the King model (Meylan et al. 2001). Considering 
these observations, the reasonable Mb of dE,N is estimated to 
be ~ -15 mag. Given a value of Mdm/Mdw, we can determine 
ttira from the above Mo-fldm relation (Burkert 1995). For conve- 
nience, fldm for Mb = -15 mag and Mdm/Md„ = 5 is referred to 
as adm,o hereafter. We mainly investigate the dE,N models with 
Mb = -15 mag, fldw = 790 pc, Mdm/Mdw = 5, adm = «dm,o, 
Mn/Md„ - 0.05, and On/fldw = 0.02. All of these values are 
reasonably consistent with observations (BinggeU & Cameron 
1991; Ferguson & Binggeli 1 994) . We also investigate the mod- 
els with fldm = 0.25adm,o to clarify the importance of the dark 
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matter halo structure of a dE,N in the formation processes of 
Gl. 

M3 1 is assumed to have the disk mass of 7.8 x 10"^ M©, and 
the bulge-to-disk-ratio of 0.25, and the maximum rotational ve- 
locity of 260 km s"', all of which are consistent with obser- 
vations (e.g., van den Bergh 2000). The initial disk plane of 
M31 is set to be the x-y plane of a simulation. The dE,N or- 
biting M3 1 has the initial position of (x, y, z) = (/?apo cos 6, 0, 
/?apo sin 6) and the initial velocity of (Vx, Vy, Vz) = (0, aVc, 0), 
where /?apo> ^> and Vc are the apocenter of the orbit, the inclina- 
tion angle with respect to the M3rs disk (i.e., the x-y plane), 
the circular velocity at the apocenter, and the parameter (0 < a 
< 1) that determines the orbital eccentricity represented by gp 
(i.e, the larger a is, the more circular the orbit is). Gl is ob- 
served to have the projected distance of 40 kpc from M3 1 and 
radial velocity of -3 1 km s~^ with respect to M3 1 (e.g., Meylan 
et al. 2001). Guided by these observations, we investigate the 
models with 6 = 30° and 80°, R^^ = 40, 80, and 160 kpc, and Cp 
= 0.62 {a = 0.5) and 0.18 {a = 0.9). We however describe the 
four representative and important models in this paper and the 
parameter values in these models (e.g., "fiducial model") are 
given in Table 1 . All the simulations have been carried out on a 
GRAPE board (Sugimoto et al. 1990) with the particle number 
of 90000. 



3. Results 

Figures 1 summarizes the morphological evolution of the enve- 
lope and the nucleus of the dE,N in the fiducial model, which 
shows typical behavior in morphological transformation from 
a dE,N into Gl. As the dE,N approaches the pericenter of its 
orbit, the strong global tidal field of M3 1 stretches the envelope 
of the dE'N along the direction of the dwarf's orbit and conse- 
quently tidally strips the stars of the envelope {T - 0.6 Gyr). 
The dark matter halo, which is more widely distributed than 
the envelope due to its larger core radius, is also efficiently re- 
moved from the dE,N during the pericenter passage. Since the 
envelope (and the dark matter halo) loses a significant fraction 
of its mass during the passage through the pericenter, the enve- 
lope becomes more susceptible to the tidal effects of M3 1 af- 
ter the pericenter passage. Therefore, each subsequent time the 
dwarf approaches the pericenter, it loses an increasingly larger 
fraction of its stellar envelope through tidal stripping. Finally, 
the envelope and the dark matter halo lose 98 % and 99 % of 
their initial masses, respectively, after four passages through 
the pericenter (7" = 4.8 Gyr). 
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Fig. 1. Upper four. Morphological evolution of the stellar envelope of the dE,N projected onto the x-z plane (i.e., edge-on view) 
for the fiducial model. For comparison, the M3 1 disk and bulge components are shown in the upper left panel and the disk size is 
indicated by a solid line in other three panels. The orbit of the dE,N with /?apo = 80 kpc and e = 0.62 is indicated by a dotted line 
in the upper left panel. The time T (in units of Gyr) indicated in the upper left corner of each frame represents the time elapsed 
since the simulation starts. Each frame is 164 kpc on a side. Lower four: Morphological evolution of the stellar envelope and the 
nucleus of the dE,N projected onto the x-y plane (i.e., face-on view) for the fiducial model. Each frame is 20 kpc on a side. 
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Fig. 2. The projected surface density profiles for the dark matter (solid line), the stellar envelope (dotted line), the nucleus (short- 
dashed line), and all these components (long-dashed fine) at T = Gyr (upper panel) and 4.8 Gyr (lower panel) in the fiducial 
model. 



The central nucleus, on the other hand, is only weakly in- 
fluenced by the tidal force as a result of its compact configu- 
ration. Because of its strongly self-gravitating nature, the nu- 
cleus loses only a small amount (~ 7%) of its mass and thus 
maintains its compact morphology during its tidal interaction 
with M31. As a result, a very compact stellar system with a 
negligible amount of dark matter is formed from the dE,N by 
T = 4.8 Gyr. The total nuclear stellar mass of the remnant 
within 5adw is ~ 1.5 X 10^ Mq, consistent with the observed 
mass of Gl. The mass-to-light ratio, M/Lg, decreases dramati- 
cally from ~ 10 to ~ 3 for r < Sajw within 4.8 Gyr. This result 
clearly explains why Gl is observed to have mass-to-Ught ra- 
tios that are much smaller than those observed for dE,Ns (~ 
10): galaxy threshing is most efficient in the outer regions of a 
dE,N where the dark matter halo dominates gravitationally. As 
is shown in Figure 2, both the surface density of the dark matter 
and that of the envelope drop by more than an order of magni- 
tude within 4.8 Gyr. If the initial central surface brightness of 
the dE,N is pg - 23 mag arcsec"^, then the final surface bright- 
ness of the dE,N at r = 3 kpc is about ps - 29.5 magarcsec"^. 
Such a faint, low surface brightness envelope will be hard to 
detect, even by existing large ground-based telescopes. 

The stripped stars from the envelope of the Gl progeni- 
tor dE,N form substructures and tidal tails in M31 halo region. 
Figure 3 shows that the tidal tail is a "rosette", reflecting the 
dE,N eccentric orbit (gp = 0.62). Most of the steUar compo- 
nents of the tail can be regarded as being located well outside 



the M31 inner halo with r < 40 kpc (~ three degrees). The 
fi-band surface brightness (jds) of the tail ranges from 28.7 
mag arcsec"^ to 30.4 mag arcsec"^ with a mean of 30.1 mag 
arcsec"^. These results imply that the tidal streams formed from 
the Gl progenitor dE,N would be hard to detect in previous ob- 
servations by Ibata et al. (2001) and Ferguson et al (2002) that 
mapped the M3 1 halo region within r ~ 40 kpc. 

Morphological properties of tails and substructures devel- 
oped during tidal destruction of a Gl progenitor dE,N depend 
not only on the initial orbit but also on the dark matter struc- 
ture of the dE,N. Figure 4 summarizes the following three de- 
pendences. Firstly, for the model with a smaller gp (= 0.18), in 
which the time scale of morphological transformation from the 
dE,N into Gl is much longer than that of the fiducial model, the 
shape of the developed tail is more like a spiral than a rosette. 
Secondly, the distribution of stars in the developed substruc- 
tures for the model with a smaller /?apo (= 40 kpc) appears to 
have a sharp boundary within which the stars can be located in a 
relatively homogeneous manner. Thirdly, the surface density of 
stars along the tail in the model with a more compact dark mat- 
ter halo (with a^m = 0.25 fldm.o) is much lower than that of the 
fiducial model owing to the smaller number of stripped stars. 
Formation of Gl from the dE,N does not occur in this model 
with = 0.25 Odm.o which suggests that a dE,N with the high 
central density of the dark matter (such as seen in NEW mod- 
els) are less likely to be a Gl progenitor dE,N (See Bekki et al. 
2003 for more discussions on this point). 
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Fig. 3. The final B-band surface brightness fis (mag arcsec"^) distribution projected onto the x-y plane for substructures/tails 
of MSl's stellar halo developed from tidal destruction of the dE,N in the fiducial model. For comparison with observations by 
Ferguson et al. (2002), the scale is given in units of degree. 



4. Discussion and conclusion 

If Gl originates from a dE,N with Mb = -15 mag, the mean 
metallicity of the M31 stellar halo consisting of stars stripped 
from the dE,N is roughly estimated as [Fe/H] = -0.96 (for 
B - V - 0.71 in the dE,N) using the observed metallicity- 
luminosity relation for dwarf galaxies, [Fe/H]=-3.43(±0.14) - 
0.157(±0.012)xMv (e.g., Cote et al. 2000). Because of its small 
mass (less than 10'' Mq) and low metallicity, the stellar halo 
formed from the tidal destruction of the Gl progenitor dE,N 
could not be the major component of the observed high density 

and metal-rich ([m/H] 0.5) MSI' stellar halo (Durrell et al. 

2000). Also, relatively metal-rich components ([Fe/H] — 0.7) 
of the stellar tail recently discovered in the M31 halo region 
(Ibata et al. 2001) might be less likely to be formed from the 
metal-poor stars stripped from Gl's progenitor dE,N. 

Recently Ferguson et al. (2002) discovered a stellar halo 
substructure located in the proximity of Gl (referred to as 
"Gl clump") and showed that total V-band magnitude and 
V-band surface brightness of the Gl clump can be estimated 



to be -12.6 mag and 28.5 mag arcsec"^, respectively, for rea- 
sonable assumptions of dust extinction. The present study has 

demonstrated that metal-poor ([Fe/H] 1) and low surface 

brightness (fiB ~ 30 mag arcsec"^) stellar halo substructures 
can be formed along the orbit of the Gl progenitor. We have 
also found that the stars in substructures close to the simu- 
lated Gl have orbital properties similar to those of the Gl. 
Accordingly, future spectroscopic observations on radial veloc- 
ities and metallicities of stars in the Gl clump will provide a 
new clue to the problem of whether the Gl clump can be tidal 
debris of the Gl progenitor dE,N. 

Our simulations suggest that if a dE,N has a higher cen- 
tral dark matter density (i.e., more compact core), it cannot be 
transformed into a giant globular cluster (Gl) because of the 
survival of its stellar envelope. Whether a dE,N can be trans- 
formed into a globular cluster by the M31 tidal field (i.e., by 
galaxy threshing) depends on whether the M31 tidal force is 
stronger than the self-gravitational force of the dark matter halo 
of the dE,N at the pericenter of the orbit of the dE,N. Therefore 
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Fig. 4. Final mass distributions projected onto the x-y plane at T = 4.8 Gyr for the fiducial model (upper left), the smaller gp one 
with gp = 0. 18 (upper right), the smaller with /?apo = 40 kpc (lower left), and the more compact dark matter model with adm 
= 0.25adm,o (lower right). Each frame is 240 kpc on a side. 



we can give the upper limit of the central dark matter density of 
the Gl progenitor dE,N, if we know its pericenter distance from 
M31. The possible central density of the dark matter halo of w 
Cen's progenitor dE,N is discussed in Bekki & Chiba (2003), 
based on the proper motion data of w Cen. 

Drinkwater et al. (2003) have recently discovered a new 
type of galactic object with effective radius of ~ 20 pc. My 
~ -12 mag, and velocity dispersion of ~ 30 km s ' in the 
Fornax Cluster. These ultra-compact dwarfs (UCDs) with My 
more than 1 mag brighter than that of Gl have been demon- 
strated to be formed by galaxy threshing in which dE,Ns with 
My < -16 mag can be transformed into UCDs owing to tidal 
stripping of the stellar envelopes of dE,Ns (Bekki et al. 2001, 
2003). We thus suggest that UCDs and giant globular clusters 
such as Gl and u) Cen can be regarded as the same class of stel- 
lar objects: The total mass or luminosity of a progenitor dE,N 
is the main difference between UCDs and Gl. However, the 
difference in the location on the My-cro relation between UCD 
and Gl cannot be explained simply by the galaxy threshing 
scenario (Bekki et al. 2003). 

The present study suggests that if Gl is not close to the peri- 
center of its orbit around M31, the stripped stars from the Gl 
progenitor dE,N can be distributed throughout the M31 outer 
halo region with Rp > 40 kpc. Therefore, future observations 



on M3 1 's stellar halo need to extend to the current limit of pre- 
vious surveys (Rp ~ 40 kpc; Ferguson et al. 2002) to reveal 
possible tidal streams and substructures formed from dE,N. 
Metallicity information of the possible streams and substruc- 
tures in M3rs outer halo is also useful to constrain the stellar 
population of the destroyed Gl progenitor dE,N. Thus, future 
deep, high-resolution, wide-area surveys of MSl's outer (Rp 
> 40 kpc) stellar halo by wide-field cameras on large ground- 
based telescopes (e.g., Suprime-Cam on Subaru) will enable 
us to determine whether Gl originates from an ancient dE,N 
orbiting M31. 

Acknowledgements. K.B. acknowledges the Large Australian 
Researcli Council (ARC). All tlie simulations described here were 
performed with GRAPE 5 systems at the National Astronomical 
Observatory in Japan. 
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